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We discuss an M-dimensional phantom elastic manifold of linear size L crushed into a small sphere 
of radius R <C L in JV-dimensional space. We investigate the low elastic energy states of 2-sheets 
(M — 2) and 3-sheets (M — 3) using analytic methods and lattice simulations. When TV > 2M the 
curvature energy is uniformly distributed in the sheet and the strain energy is negligible. But when 
N = M + 1 and M > 1, both energies appear to be condensed into a network of narrow M — 1 
dimensional ridges. The ridges appear straight over distances comparable to the confining radius R. 
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It has long been known that a thin elastic plate will 
develop narrow ridges under a variety of compressive 
boundary conditions. These ridges may be seen daily 
in the way a pillowcase or trouser leg deforms to wrap 
its contents, often exhibiting a diamond pattern famil- 
iar from compression studies of thin metal cylinders . 
The linear scars in a crumpled sheet of paper are also a 
record of this mechanism |2[] . Recently it was discovered 
that the structure of these ridges could be accounted for 
using linear elasticity theory, valid in the limit that the 
ridge length X is much greater than the plate thickness 
h HH . Witten and Li used a scaling argument to predict 
that the ridge width w ~ /i 1 / 3 ^ 2 / 3 and the total elastic 
energy E ~ Yh 3 (X/h) 1 ' 3 where Y is the Young's modu- 
lus H . Lobkovsky et al. verified these scaling laws using 
both numerical simulations and an asymptotic analysis of 
the Von Karman equations for a thin plate Q. His sim- 
ulations showed that the material strains and curvatures 
decay rapidly to zero in the direction transverse to the 
ridge. The length scale of this decay is the ridge width w, 
which goes to zero with the thickness of the plate. How- 
ever, these ridges were analyzed only in idealized, sym- 
metrical deformations. Their applicability to stochasti- 
cally crumpled sheets has not been explicitly shown. 

Ridge formation is a mode of spontaneous condensa- 
tion of energy into a small subset of the available vol- 
ume. As such it resembles the spontaneous organization 
of dislocations into grain boundaries in a strained crystal 
or the formation of Prandtl boundary layers in laminar 
fluid flow §,§. 

To understand the necessary conditions for this con- 
densation, and its consequences, it is useful to consider 
the general problem of an M-dimensional elastic mani- 
fold crushed by a hypersphere in iV-dimensional space. 
In Ref . @ we derived the elastic energy functional of the 
manifold by considering the small thickness limit of an 
iV-dimensional elastic solid with an extent 0(L) in M di- 
rections and a thickness h in N—M transverse directions. 
Using this formalism we found that a deformed hypersur- 
face (M = N — 1) under a specific boundary condition 
exhibits a ridge with scaling properties analogous to those 



of Ref. H . In this Letter we explore the more general sit- 
uation of an M-sheet with a free boundary confined by a 
small sphere. Using analytic and numerical methods we 
demonstrate two complementary behaviors: 

Case 1 (N > 2AI ) The energy E c associated with cur- 
vature is distributed uniformly in the manifold and the 
energy E s associated with strain is negligible by compar- 



Case 2 (N = M + 1 > 2) The strain energy and curva- 
ture energy are comparable, with E s jE\, ~ 0.2. Both en- 
ergies are condensed into a network of narrow, M — 1 
dimensional ridges as described in Refs. j2l and Qj. 

In the remaining cases, namely M + 1 < N < 2M, we 
anticipate that the strain energy will also be comparable 
to the curvature energy ||. We have not investigated 
these cases in detail. 

The differential geometry of the deformed manifold 
is conveniently discussed using a Euclidean coordinate 
patch on the flat manifold (x Q |a S [l,Af]), called the 
manifold coordinate system ||. Any deformation may 
then be written as a map r{x) from the manifold coor- 
dinates to A^-dimcnsional Euclidian space. The strain 
tensor is u a p(x) = (l/2)(d a r ■ dpr — 8 a p) and the ex- 
trinsic curvature tensor is K a p — d a dpr ||. Under the 
usual assumptions of linear elasticity theory (u a p -C 1 
and <9 7 u Q/ 3 -C \Kn V \ |1C|] ) the elastic energy func- 

tional becomes 



E = J dx M {(i ((u al3 ) 2 + c (u aa ) 2 ) 



K I K a p ■ K a p + c K aa ■ K t 



»)} 



(1) 



where /j, ~ Yh N ~ M is an effective Lame coefficient, 
k ~ Yh N ~ M+2 is an effective bending rigidity, and cq 
is a dimensionless constant [Q,|llJ. Summation over re- 
peated indices is implied. We shall refer to the term 
quadratic in the strains as the strain energy E s , and the 
term quadratic in the curvatures as the curvature energy 



1 



E c . We consider only phantom (not self- avoiding) mani- 
folds in this paper, since we don't expect self-avoidance 
to alter the qualitative conclusions. Our compressive 
boundary condition is a frictionless, hard-walled sphere 
with an initial radius Ri > L. To crush the manifold the 
radius is slowly decreased to a final value R -C L. 

In our simulations we model a thin clastic 2-sheet 
(M = 2) using a triangular network of nodes connected 
by springs (see Fig. 1), following the work of Seung and 
Nelson . The discretized strain energy is the sum 
over the Hooke's law energy of each spring 

E s =Y J \ c s{h~h) 2 (2) 

i 

where c s is a spring constant controlling the strain energy 
U is the length of spring i, and Iq is the lattice constant. 




b 

FIG. 1. A portion of the triangular lattice and the 
simple cubic lattice used to discretize a 2-sheet (M = 2) 
and a 3-sheet (M = 3) respectively. 



To discretize the curvature energy it is convenient to 
begin by defining the normal tensor for a two-dimensional 
lattice in iV-dimensional Euclidean space. If a(x) and 
b(x) are two independent vectors tangent to the 2-sheet 
at x, then we may define a normal tensor via 

(n[a,b])i 1 ...i N _ 2 — ^i 1 ---i N 0'i N ^ 1 bi N , (3) 

where and bi are the Cartesian components of a and 
b and e is the antisymmetric Levi-Civita tensor. In 3- 
space this reduces to the usual normal vector n = a x b. 
Defining the inner product 

n i ' n 2 = (jv _ j^nj ( n i)n---»«-M( n 2K---i w - M (4) 

with M = 2, the unit normal tensor may be expressed 
n = n/v^n • n. By substituting a — d\v and b — d^r, it 
is straightforward to prove that d a n ■ d a n = K a p • K a p. 
The curvature energy can therefore be represented as the 
square of the discrete derivative 



Ec = J2 \cb\Hi - ft/ (5) 

<y> 

where Cf, is a bending constant, fij is the unit normal of 
triangle i (see FIG. 1), and the sum is taken over each 
pair of triangles which share a common edge. Compari- 
son with Eq. (|l|) gives /i ~ c s /Iq, k ~ Q,/Zq, and Co = 
fL2f . The effective thickness of the 2-sheet is determined 
by the ratio Cb/c s ~ h? . We take (c^/cs) 1 / 2 > l /6 to 
minimize lattice effects in our simulations. 

The discrete model of an elastic 3-sheet (M = 3) is 
qualitatively similar. The 3-sheet is a simple cubic lat- 
tice of nodes connected by springs. As shown in Fig. 
1, there are edge springs with constant c e and length Iq 
and diagonal springs with constant Cd and length \/21q. 
The requirement that the lattice be elastically isotropic 
for small strains is q = 2c e . To calculate the curvature 
energy we first define the normal tensor on a 3-manifold 
in TV-space 

(n[a,b,c]) il ... lN _ 3 = e il ...i N a iN _ 2 bi N _ 1 Ci N (6) 

where a{x), b(x), and c(x) span the tangent space of the 
manifold at (x). Using the inner product Eq. (^) with 
M = 3, we again have d a n ■ d a n — K a p ■ K a p. Each 
unit cube is conceived as five adjacent tetrahedra with 
the springs as edges. Each cube then has four "corner" 
tetrahedra surrounding one "middle" tetrahedron. The 
sum in Eq. (^|) is taken over all pairs of tetrahedra that 
share a common triangular face (see Fig. 1). Lastly, we 
relate the bending constant c cc for a pair of corner tetra- 
hedra to the bending constant c cm for a corner-middle 
pair. Isotropy requires c cm = 2c cc . 

The initial condition in our simulations is a hexagonal 

2- sheet (M = 2) or a roughly spherical 3-sheet (M = 3) 
with a small random displacement added to the positions 
of the nodes We model a hyperspherical container 
of radius tq with the potential V sp here — Ei( r i/ r o) 12 ) 
where is the distance from the origin to node i. We 
use the rms radius of the manifold as our measure of the 
confining radius R: R 2 — J2i r ? I (J2i !)• To crush the 
manifold, the value of r is repeatedly decreased, and 
the energy of the manifold minimized at each step using 
a conjugate gradient routine E3|. Using this method a 
hexagonal 2-sheet with long diameter L — 16(Mo may be 
crushed to a radius r — 20l in three days of CPU time 
on an IBM RISC 6000. An approximately spherical ball 
with diameter L = 301q may be crushed to ro = 6.0/o in 
about five days. 

The simplest instance of Case 1 is a thin elastic 
rod crushed within a circle in 2-space. As the radius 
of the circle is decreased the rod buckles and devel- 
ops a curvature of 0(1 /R). Because the rod is free 
to reptate parallel to its length, its strain energy re- 
mains identically zero. Analogously, a thin plate in 4- 
space can curl without strain into an arbitrarily small 

3- sphere. One possible embedding which demonstrates 
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this is r = p[cos(x\/ p), sm(x\/p), cos(x2/p), sin^/p)] 
where p = R/y2. Note that r 2 = R 2 everywhere. One 
can verify by substitution that the strain tensor u a p 
is identically zero, |J?n(x)| = \K 2 2{x)\ = \/2/R, and 
Etot — k(L/R) 2 . The key feature of this embedding is 
the separation of the manifold coordinates into indepen- 
dent, 2-dimensional subspaces of 4-space. Analogously, 
whenever N > 2M we may write the embedding 

r = p[cos(xi/p),sm(xi/p),cos(x 2 /p),sm(x 2 /p), 

. . . ,cos(x M /p),sm(x M /p),0, 0, . . .] (7) 



a network of narrow ridges which connect the vertices. 
The strain energy density is similarly localized. The ra- 
tio E s /E c w 10%. This is consistent with the devia- 
tions from the asymptotic value E s /E c = 1/5 reported 
by Lobkovsky for short ridges X w lOh 0]. 



where p — R/vM. This embedding is an existence proof 
that the manifold need not strain during the compres- 
sion. It is also the global minimum of the elastic energy 
at fixed R (neglecting corrections within R of the man- 
ifold boundary). This deformation removes the isotropy 
of the manifold whenever M > 1. Thus, e. g. in order to 
return to the starting position rone must move y/2 as far 
in the (1, 1) direction as in the (1, 0) direction. With this 
broken symmetry comes a degeneracy. For every distinct 
orientation of the coordinate system {x\, Xm}, Eq. ([?]) 
yields a new minimum-energy configuration. 

Our simulations of a 2-sheet in 4-space and a 3-sheet 
in 6-space behave as anticipated for Case 1. In both in- 
stances we observe the minimum energy embedding dis- 
cussed above, with a curvature energy density uniform to 
10% and negligible strain energy E s < 10 _3 i? c . The only 
significant deviation from the ideal embedding is that the 
manifold flattens out within R of the manifold edge. As 
a result, the curvature energy density decays to zero near 
the edge and there is a systematic downward correction to 
the total energy of the form E tot ~ [L— (const. )R) M /R 2 . 
Depending upon the initial condition, we also observe 
several metastable states with a nonuniform curvature 
energy density and a higher total energy. As compression 
proceeds these metastable states make the transition to 
the minimum energy embedding. 

Our understanding of Case 2 deformations is incom- 
plete. It is based on previous analytic and numerical work 
for a variety of simple deformations, including some on 
general hypersurfaces plfjj . These cases all show the con- 
densation of the elastic energy into a fraction O(hfX) 1 ^ 3 
of the manifold volume, suggesting that a generic com- 
pression yields similar condensation. The simple cases 
also exhibit a strict proportionality between strain and 
curvature energy: \im^ h / X )-^a E s = E c /5. 

We simulated Case 2 deformations using a 2-sheet in 
3-space and a 3-sheet in 4-space. For the 2-sheet we see 
the anticipated spontaneous formation of a network of 
linear ridges. Fig. 2 shows the curvature energy density 
in the manifold coordinates of a hexagonal 2-sheet with 
long diameter L — 160^o crushed to a radius R = 271q. 
About 40% of the total curvature energy occupies just 
2% of the total area. This is due to the formation of 
point vertices in the 2-sheet. These vertices are the tips 
of cone-like deformations. The next 40% occupies 20% 
of the area. We see that this energy is condensed into 




FIG. 2. Curvature energy distribution in a hexagonal 
2-sheet with moduli c s = 1.0 and c& = .025 and long di- 
ameter L = 160^o crushed to a radius tq — 27Zq. Darker 
regions have higher energy density. The crumpled con- 
figuration is shown in the lower right. 

Fig. 2 suggests that the ridges forming the network 
are of comparable length. This length is similar to the 
diameter of the confining sphere shown in the lower right 
corner. If we make the scaling hypothesis that the mean 
ridge length X ~ R, then the number of ridges in the 
2-sheet scales like (L/R) 2 . Since the energy of a single 
ridge scales as Yh 3 {R/h) 1 ^ 3 JO, the total elastic energy 
E tot should obey E tot ^ YhLr(h/R) 5 / 3 . Our simulation 
results are consistent with E to t ~ R~ 5 ^ 3 However, 
due to boundary effects and the small va ues of the as- 
pect ratio L/h< 300 the evidence for this scaling is only 
suggestive. 

Our simulations of an elastic 3-sheet in 4-space at- 
tained aspect ratios of only L/h ~ 30. Since the ridge 
properties derived in Ref. (?]] are well-defined only in the 
limit that the ridge is much longer than h, these simu- 
lations provide only qualitative support for Case 2. The 
simulations show the anticipated concentration of strains 
and curvatures into linear and planar structures. Fig. 
3 shows two curvature energy isosurfaces in the man- 
ifold coordinates of a 3-sheet with diameter L = 30lo 
crushed to r = 12Z - The upper isosurface encloses 65% 
of the total curvature energy and 16% of the volume. We 
see distinct planar structures. The lower isosurface en- 
closes 23% of the curvature energy and only 3% of the 
volume. The energy is concentrated on the set of lines 
where planes intersect. The lines are the analogs of the 
point-like vertices in FIG. 2. The ratio E s /E c is 15%. 

We anticipate scaling behavior for the energy of a 
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crushed hypersurface analogous to that for a crushed 2- 
sheet. If the linear size of the ridges is roughly the di- 
ameter of the confining sphere, we may generalize the 
argument for the 2-sheet to infer E tot ~ YhL M (h/R) 5 / 3 
. The limitations of our simulation prevented us from 
testing this prediction. 



ganizations of matter - an elastic manifold. Increasingly, 
the behavior of simple manifolds in general spatial dimen- 
sions is invoked to account for fundamental processes |jl7) . 
The symmetry-breaking and ridge-forming mechanisms 
explored here may prove relevant for understanding such 
behavior. 




FIG. 3. Curvature energy isosurfaces in the manifold 
coordinates of an elastic 3-shcct (L = 30lo) in 4-space 
crushed to a radius ro = 131q. The top (bottom) surface 
encloses 65% (23%) of the total energy in 16% (3%) of 
the volume. 

This work has implications for real crumpled sheets. 
It is the first exploration of the distribution of energy 
in these sheets and it offers qualitative support for the 
model of a crumpled sheet as a network of ridges . The 
possibility of crushing without strain in high dimensions 
may be relevant for studies of thermal crumpling in gen- 
eral dimensions 16 1. More broadly, this work identifies a 
new mechanism of energy condensation into a small sub- 
space of an available volume. This condensation happens 
in arbitrary spatial dimensions in one of the simplest or- 
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